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Antiferromagnetism of Ni2NbB06 with S = 1 dimer quasi-lD armchair chains 
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Long range antiferromagnetic (AFM) ordering of Ni spins in Ni2NbB06 has been studied with 
single crystal from spin susceptibility measurement and compared with the ah initio calculation 
results consistently. Below Tn~ 23.5 K, the *9 = 1 spins align along the a-direction for edge-shared 
NiOe octahedra which form crystallographic armchair chains along the ^-direction. The isothermal 
magnetization M(H) below Tn shows spin-flop transition for magnetic field above ~ 36 kOe along 
the a-axis, which indicates the spin anisotropy is along the a-direction. The electronic and magnetic 
structures of Ni 2 NbB 06 have also been explored theoretically using density functional theory with 
generalized gradient approximation plus on-site Coulomb interaction (U). These calculations sup¬ 
port the experimentally observed antiferromagnetism of Ni2NbB06. In particular, the long range 
AFM ordering below Tn can be dissected into armchair chains which consists of S' = 1 dimers of 
J2~2.43 meV with ferromagnetic (FM) intra-chain and inter-chain couplings of size J 2 I. 

PACS numbers: 75.50.Ee, 75.30.Et, 75.30.Gw, 71.15.Mb 


I. INTRODUCTION 

Low dimensional magnetic materials have attracted 
considerable attention due to their interesting low tem¬ 
perature properties with the involved strong quantum 
fluctuations^. Extensive studies of materials with ge¬ 
ometric frustration on square, triangular, zigzag chains 
and zigzag ladders spin systems have been explored for 
the diverse magnetic ground states. The zigzag spin 
chain of S = 1/2 with antiferromagnetic (AFM) inter¬ 
actions between nearest neighbor (NN) and next nearest 
neighbor (NNN) is about the most commonly studied 
frustrated systenP^. In zigzag spin chain system with 
5 = 1, the ground state phase diagram as a function of 
anisotropy and ratio between NN and NNN interactions 
exhibits different phase^^. Metal borates are expected 
to be good candidates to serve as links for transition 
metal pol^edra giving rise to different low-dimensional 
structured. Another important role of the borate anions, 
being non-magnetic, is to allow transmission of magnetic 
interactions via a super-superexchange routd^^^. 

In the present work, we report the crystal growth and 
the magnetization measurement results along the three 
principal directions of Ni2NbB06. Grystallographically 
Ni2NbB06 has been found to be a 5 = 1 armchair spin 
chain systerrP^. We found that a long range AFM spin 
ordering exists below Tn^ 23.5 K. A sizable inter-chain 
coupling leads to the 3D long range AFM spin ordering 
with an on-site anisotropy along the a-direction, which 
is as confirmed by the field-induced spin flop transition. 
We also studied the electronic and magnetic properties of 
Ni2NbB06 within the density functional theory with the 
generalized gradient approximation. We found that the 
system consists of unconventional armchair chains which 


are formed with ferromagnetically coupled 5 = 1 dimers 
with intra- and inter-chain coupling constants which are 
nearly half of that for the 5 = 1 dimer. An interpreta¬ 
tion on the experimental observation on the AFM and 
spin flop transition is provided and compared with the 
calculated results. 


II. EXPERIMENTAL AND COMPUTATIONAL 
DETAILS 

Single crystal of Ni2NbB06 was grown by a flux 
method using borate as the solvent. A mixture of 6.6 g 
NiO, 20 g Nb205 and 33 g of Na2B407 were placed in 
a platinum crucible and heated to 1250°G in a box fur¬ 
nace for 24 hours. The furnace was slowly cooled down 
to 850°C at a rate of 3°C/h and then cooled down to 
room temperature at the rate of 80°C/h. The single 
crystals in green color (shown in the inset of Fig. [^b)) 
were separated from the borate flux by leaching with 
a dilute solution of HNO3. The crystal structure and 
phase purity of the samples were checked by powder X- 
ray diffraction (XRD) using the synchrotron X-ray of 
A = 0.619 A(NSRRC, Taiwan) at room temperature. The 
field cooled (FC) and zero field cooled (ZFG) magneti¬ 
zation curves were measured in a commercial Vibrating 
Sample Magnetometer (VSM, Quantum Design, USA) 
from 1.8 K to 300 K in the presence of various applied 
magnetic fields. The isothermal magnetization (M) data 
were also recorded at selected temperatures. 

Theoretical calculations have been performed based on 
first-principle density functional theory (DFT) with gen¬ 
eralized gradient approximation (GGA)^. The on-site 
Coulomb energy U has been taken into account using 
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FIG. 1. (Color online) (a) The powder XRD pattern from the 
pulverized as grown Ni2NbB06 crystal, (b) The XRD pattern 
of Ni 2 NbB 06 crystal perpendicular to the large surface to 
show (hOO) peaks with preferred orientation. Inset is the as 
grown single crystal. 


the GGA+U schem^i^. We have used effective Ueff = 
{U — J) = 6 eV for the Ni atoms in the GGA+U cal¬ 
culations. We used the accurate full-potential projector- 
augmented wave (PAW) methocP^ implemented in the 
Vienna ab initio simulation package (VASP)^^^^. Ex¬ 
perimental lattice parameters were used in the calcula¬ 
tion. The primitive unit cell contains four Ni2NbB06 
formula units. In the present calculations, we used the 
tetrahedron method with Blochl corrections for the Bril- 
louin zone integration with a F-centered Monkhorst-Pack 
k-point mesh of (8 x 10 x 18). A large plane wave cutoff 
energy of 500 eV was taken, and the convergence criterion 
for the total energy was 10“^ eV. 

III. RESULTS AND DISCUSSION 
A. Crystal structure 

The powder XRD pattern of the polycrystalline sample 
obtained from the pulverized as-grown Ni 2 NbB 06 single 


crystal sample is shown in Fig.j^a). All diffraction peaks 
can be indexed to the orthorhombic structure with space 
group Pnma, without any observable trace of impurity 
phase. The structural parameters were refined by the 
Rietveld technique with good quality refinement param¬ 
eters (Rwp = 1.67 % and R^ = 1.04 %). The obtained 
values of the lattice parameters are a = 10.0690(1) A, 
b = 8.6266(2) A, and c = 4.4932(3) A, which are in ^od 
agreement with previously reported valueFig. mb) 
illustrates the single crystal XRD pattern with peaks in¬ 
dexed for the preferred orientation perpendicular to the 
(hOO) planes. This compound could also be viewed as lay¬ 
ers containing armchair chains of edge-shared NiOe octa- 
hedra, where each pair of NiOe along the 5-direction are 
edge-shared with both NbOe octahedra and BO 4 tetra- 
hedra, as illustrated in Fig. 



FIG. 2. (Golor online) (a) A layer of edge-shared NiOe octahe¬ 
dra (green) shown in c-projection containing armchair chains 
along the 5-direction, where every BO4 tetrahedron (blue) in 
the neighboring layer bridges a pair of NiOe and one NbOe 
octahedra (yellow) through face-sharing, as shown in (b). 


B. Magnetic susceptibility 

Fig. shows the temperature dependence of magnetic 
susceptibility x( T) and the corresponding inverse suscep¬ 
tibility x“^(r) measured at an applied magnetic field of 
10 kOe in the temperature range 2-300 K for pulverized 
powder of the as-grown Ni2NbB06 crystal. The x(T) 
curve shows a Curie-Weiss like behavior at high temper¬ 
ature and a sharp peak is observed at 24 K, indicating 
the onset of a antiferromagnetic ordering. The ordering 
temperature Tn = 23.5 K is defined by the sharp peak 
through d(x7")/dr. At T > 50 K, the x(^) data can be 
fitted with the Curie-Weiss law (x(^) = Cy(r-6>)) satis¬ 
factorily using the Curie constant C= 1.31 and the Curie- 
Weiss temperature {0 = 9.5 K), as shown by the red solid 
line in Fig. The effective moment of /ie// = 3.23 
per Ni^+ extracted from Curie constant is higher than the 
expected spin-only value of jiicaic = 2.83 fiB for S = 1, 
which suggests the existence of a partially unquenched 
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FIG. 3. (Color online) The temperature dependence of mag¬ 
netic susceptibility x(^) ^md the corresponding inverse sus¬ 
ceptibility x“^(T) measured at an applied magnetic field of 
10 kOe in the temperature range 2 — 300 K for the pulverized 
as-grown Ni 2 NbB 06 crystals. 


orbital contribution. The fitted value of ^ = 9.5 K sug¬ 
gests the existence of an average ferromagnetic (FM) cou¬ 
pling among spins at high temperature but AFM order¬ 
ing occurs at T7v^23.5 K, which suggests that the mag¬ 
netic interactions must consider couplings beyond nearest 
neighbor spins and have different signs containing both 
ferromagnetic (FM) and AFM couplings, as verified later 
by our ab initio studies in the following. The x(T) data 
above 200 K were also fitted by high temperature series 
(HTS) expansion up to 8^^ order.^^ The fitting parame¬ 
ters are found to be ^ = 2.02 and the exchange interaction 

(J/k^) = -6.6 K. 
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FIG. 5. (Color online) (a) Isothermal magnetization M(H) 
curves of Ni 2 NbB 06 for the magnetic field parallel to differ¬ 
ent crystallographic axis (b) dM/dif curves of Ni 2 NbB 06 for 
the magnetic field parallel to a-axis at some selected temper¬ 
atures. 



T(K) 


FIG. 4. (Color online) The magnetic susceptibilities 
x(T) for Ni 2 NbB 06 single crystal measured in an applied 
magnetic field of 10 kOe parallel to all three crystallo¬ 
graphic axis. The average of single crystal susceptibility 
(x( T) = (Xa+X6+Xc)/3) agrees perfectly with those measured 
using powder sample directly. 

Anisotropic magnetic susceptibilities x(T') for 


Ni 2 NbB 06 single crystals were measured in an applied 
magnetic field of 10 kOe parallel to all three crystallo¬ 
graphic axes, as shown in Fig. There is no deviation 
between Xzfc{T) and Xfc{T) throughout the measured 
temperature range. Below Tn, the anisotropy becomes 
significantly enhanced as shown in the inset of Fig. 
which indicates that the spins are aligned along the 
a-axis for in the 3D AFM long range ordering. 

To understand the AFM behavior of Ni2NbB06 bet¬ 
ter, we measured magnetization as a function of magnetic 
field H along the two crystal orientations above and be¬ 
low Tn, as shown in Fig.|^a), where no field or tempera¬ 
ture hysteresis were observed. For magnetic field parallel 
to the a-axis, the magnetization reveals jump at a critical 
field He near Hq = 36.7 kOe below T^v at 2 K, which 
is attributed to a spin-flop transition when the spin sus¬ 
ceptibility changes abruptly to a higher level, z.e., the or¬ 
dered spins originally aligned along the a-direction flop to 
the direction perpendicular to the external field. As ex¬ 
pected for the ordered spins aligned along the a-direction, 
we do not observe spin-flop transition manifested in M(H) 
for H is applied parallel to b- or c- axis. We may summa¬ 
rize the H-T phase diagram for Ni2NbB06 based on the 
magnetic field and temperature dependence of M{T,H) 
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FIG. 6. (Color online) H - T phase diagram of Ni2NbB06 for 
magnetic field parallel to a-axis obtained from magnetization 
measurements. The shaded area is related to the spin-flop 
(SP) transition defined by the start and end points of the 
transition from Open circles represent measured Tn 

and solid line is guide to the eye. 

with field applied parallel to a-axis, as shown in Fig. 

A small field dependence of Tn is also shown, where the 
boundary of spin-flop transition is indicated according to 
the onsets of dM/diF peaks shown in Fig. [^b). 


C. Theoretical calculations 

Within first-principle density functional theory, we 
first calculated the total energy (Efm) for the ferromag¬ 
netic state. The total energy per formula unit (f.u.) is 
-72.7578 eV. In order to find out the magnetic ground 
state of the system, we have considered various magnetic 
configurations possible within the unit cell. Three config¬ 
urations corresponding to the possible magnetic ground 
states have been used to estimate the exchange interac¬ 
tions, and three coupling constants (J^) are considered 
based on the three shortest Ni-Ni distances, as shown in 
Fig.0 Here Ji, J 2 , and J 3 represent the exchange cou¬ 
plings between two neighboring Ni atoms corresponding 
to the Ni-Ni distances of 2.987 A, 3.099 A, and 3.436 A, 
respectively, as shown in Fig. Both ferromagnetic as 
well as antiferromagnetic alignments of Ni spins are con¬ 
sidered and labelled as configurations A, B, and C. The 
calculated total energies of all these configurations are 
summarized in Table [J . 

We find that the configuration A (Fig.[^ has the lowest 
energy, therefore configuration A is the expected mag¬ 
netic ground state of the system. In this configuration, 
all the NN Ni ions within the armchair chain along the 
6 -direction are antiferromagnetically coupled as a 5 = 1 
dimer, and the other two NNN Ni ions, z.e., the inter¬ 
dimer within the armchair chain and the inter-chain con- 




FIG. 7. (Golor online) A schematic representation of the spin 
arrangement for configuration A. 


plings are ferromagnetically coupled. 

The calculated magnetic moment ofNi is ^1.78 fiBi 
which is slightly smaller than the expected value of 2 /i^ 
for Ni^+, suggesting that some of the magnetic moment 
lies outside the nickel atomic sphere used. Fig. shows 
the band structure (top panel) and density of states (bot¬ 
tom panel) of configuration A. We find that the system 
has a large gap >3 eV, which indicates that the system 
is insulating and consistent with the experimental obser¬ 
vation. It is quite clear from the site-resolved density of 
states that the valence band is mainly composed of nickel 
3d and oxygen 2p states. Thus the magnetic structure 
should be primarily decided by the spin-exchange cou¬ 
pling via Ni-O-Ni. 

Next we evaluated the intra-chain and inter-chain mag¬ 
netic exchange couplings among Ni spins. The exchange 
interaction between all the nearest-neighbor Ni ions along 
the 6 -direction in the a 6 -plane is denoted as J 2 , and the 
other two couplings are Ji (intra-chain) and J 3 (inter¬ 
chain), as illustrated in Fig. In order to evaluate the 
exchange couplings, we have considered the obtained to¬ 
tal energy of the unit cell of Ni 2 NbB 06 as the sum of the 
nearest-neighbor spin-spin interactions in terms of the 
spin Heisenberg model H = Eq JijCfi • crj. Here 

Jij is the exchange interaction parameter between the 
nearest-neighbor Ni site i and site j, and a[ (crj) is the 
unit vector representing the direction of the local mag¬ 
netic moment at site i (j). The total energy per unit cell 


TABLE I. Galculated total energy AE (relative to the total 
energy of FM state Efm = 72.7578 eV/f.u.), total magnetic 
moment atomic moment of Ni and band gap Eg 


Gonfig. 

AE 

(meV/f.u.) 

(Ats/f.u.) 

(/is/atom) 

Eg 

(eV) 

FM 

0.0 

4.0 

1.78 

2.9 

A 

-4.87 

0.0 

1.78 

3.3 

B 

-2.0 

0.0 

1.78 

3.3 

G 

-2.16 

0.0 

1.78 

3.3 
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FIG. 8. (Color online) Band structure (top panel) and density 
of states (bottom panel) of configuration A. Top of the valence 
band has been set to zero. 


TABLE II. Calculated exchange interaction parameters (Ji) 
and the corresponding nearest-neighbor Ni-Ni distances. 


Ji J2 J3 
T(meV) -1.43 2.43 -1.27 
Ni-Ni(A) 2.987 3.099 3.436 


for all considered magnetic configurations are given by 


Efm = ^o+4(J3 +Ji)+4J2, Ea = F^o+4(J3 +Ji) — 4 J 2 , 
Eb = F^o+4 (J3 -Ji)-4 J2, and Ec = EQ-AJ 2 . Solving 
the above mentioned equations we get the values of all 
exchange interactions listed in Table |n| , where J > 0 for 
AFM interaction and J < 0 for FM interaction, and the 
constant Eq contains all spin-independent interactions. 
Although the Ni-Ni distances for all these FM and AFM 
couplings are very close, the strength of the AFM cou¬ 
pling is nearly twice of the FM coupling, which could 
be due to the strong bonding of boron tetrahedra that 
bridge the NiOe pairs within each armchair chain. The 
AFM superexchange coupling has the largest magnitude 
of 2.43 meV, which roughly corresponds to rN^28 K and 
in good agreement with the experimental observation of 
Tn-23.5 K. 

IV. SUMMARY 

The magnetic properties of Ni2NbB06 containing arm¬ 
chair chains have been studied in detail through M(H,T) 
measurement of the single crystal sample and com¬ 
pare with models predicted by the ab initio calcula¬ 
tions. A long range AFM spin ordering observed be¬ 
low rN^23.5 K. The spin-flop transition of critical field 
36.7kOe at 2 K is found along a-axis and the H-T 
phase diagram is constructed accordingly. Within first- 
principle density functional theory, we have calculated 
the electronic and magnetic structures with exchange in¬ 
teractions that agree satisfactorily with the experimen¬ 
tal results. We have established that Ni2NbB06 con¬ 
sists of unusual armchair chains which are formed with 
5 = 1 dimers with ferromagnetic intra- and inter-chain 
couplings. 


V. ACKNOWLEDGEMENT 

FCC acknowledges the support provided by MOST- 
Taiwan under project number MOST 102-2119-M-002- 
004. GYG acknowledges the financial support for this 
work from the Academia Sinica Thematic Research Pro¬ 
gram and the Ministry of Science and Technology of Tai¬ 
wan. 


* gyguo@phys.ntu.edu.tw 

^ fcchou@ntu.edu.tw 

^ H. T. Diep, ’’Magnetic systems with competing interac¬ 
tions (Frustrated spin systems)”. World Scientihc, Singa¬ 
pore, 1994. 

^ M. Troyer, H. Tsunetsugu and D.Wurtz, Phys. Rev. B 50, 
13515 (1994). 

^ S. R. White and I. Affleck, Phys. Rev. B 54, 9862 (1996). 

^ M. Matsuda and K. Katsumata, J. Magn. Magn. Mater. 
140-144, 1671 (1995). 

^ T. Hikihara, M. Kaburagi, H. Kawamura and T. Tone- 
gawa, J. Phys. Soc. Japan 69, 259 (2000). 


® A. Kolezhuk, R. Roth and U. Schollwock, Phys. Rev. Lett. 
77, 5142 (1996). 

^ J. D. Grice, P. G. Burns and F. G. Hawthorne, Gan. Min¬ 
eral. 37, 731 (1999). 

® V. V. Grishachev, O. S. Kolotov, A. P. Krasnojon and V. 

A. Pogozhev, J. Magn. Magn. Mater. 241, 81 (2002). 

® J. P. Attheld, A. M. T. Bell, L. M. Rodriguez-Martinez, 
J.M.Greneche, R. J. Gernik, J. F. Glarke and D. A. 
Perkins, Nature 396, 655 (1998). 

J. G. Fernandes, F. S. Sarrat, R. B. Guimaraes, R. S. Fre¬ 
itas, M. A. Gontinentino, A. G. Doriguetto, Y. P. Mas- 
carenhas, J. Ellena, E. E. Gastellano, J-L. Tholence, J. 
Dumas, and L. Ghivelder, Phys. Rev. B 67, 104413 (2003). 































6 


J. C. Fernandes, R. B. Guimaraes, M. A. Continentino, R. 
Rapp, J-L. Tholence, J. Dumas, Y. Blancquaert, S. Yates, 
and C. Paulsen, Phys. Rev. B 69, 054418 (2004). 

G. B. Ansell, M. E. Leonowicz, M. A. Modrick, B. M. 
Wanklyn and F. R. Wondre, Acta Gryst. B 38, 892 (1982). 
J. P. Perdew, K. Burke and M. Ernzerhof, Phys. Rev. Lett. 
77, 3865 (1996). 

S. L. Dudarev, G. A. Botton, S. Y. Savrasov, G. J. 
Humphreys and A. P. Sutton, Phys. Rev. B 57, 1505 
(1998). 


P. E. Blochl, Phys. Rev. B 50, 17953 (1994). 

G. Kresse and J. Hafner, Phys. Rev. B 47, 558 (1993); 49, 
14251 (1994). 

G. Kresse and J. Furthmiiller, Gomput. Mater. Sci. 6, 15 
(1996). 

G. Kresse and D. Joubert, Phys. Rev. B 59, 1758 (1999). 
G. S. Rushbooke and P. J. Wood, J. Molec. Phys. 1, 257 
(1958). 



